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The neutrino driven wind during a core collapse supernova is an attractive site for r-process 
nucleosynthesis. The initial electron fraction Ye in the wind depends on observable neutrino energies 
and luminosities. The mean antineutrino energy is limited by supernova SN1987A data while lepton 
number conservation constrains the ratio of antineutrino to neutrino luminosities. If Ye, in the wind, 
is to be suitable for rapid neutron capture nucleosynthesis, then the mean electron neutrino energy 
could be significantly lower then that predicted in present supernova simulations. Alternatively, 
there could be a rapid increase of the mean antineutrino energy at late times. However, this is not 
seen in the SN1987A data. Finally, there could be new neutrino physics such as oscillations to sterile 
neutrinos. 



The neutrino driven wind above a protoneutron star in a core collapse supernova is an attractive site for r-process 
nucleosynthesis. In the r-process, seed nuclei rapidly capture free neutrons to produce about half of the heavy elements 
10,^. Many simulations, for example have explored a range of physical conditions including entropy, expansion 
time scale, and electron fraction (number of electrons or protons per baryon) Yg necessary to produce r-process 
elements with solar system abundances. 

The inital electron fraction Yg is an important parameter that determines the number of free neutrons. If there 
are too few free neutrons per seed nucleus, then the heaviest elements may not be produced. A reasonable minimum 
requirement for an r-process, producing solar system like abundances, is that Ye be less then 1/2. If Ye is greater then 
1/2, all of the neutrons may be quickly incorporated into alpha particles leaving only free protons. 

In a core collapse supernova, high neutrino luminosities eject some baryons from the surface of the protoneutron 
star into a neutrino driven wind. Many authors have explored r-process nucleosynthesis in this wind |^-^. A number 
of detailed recent simulations have had difficulties reproducing a successful r-process in this wind, see below. In this 
paper we express some of these difficulties in a more model independent form by using data from Supernova Snl987A 
and plausible physical assumptions to constrain the r-process. It is hoped this model independence will help clarify 
assumptions and problems in the search for the r-process site. 

The electron fraction Ye in the wind is set by the relative rates of the neutrino capture reactions, 

i^e + n^p + e', (1) 



i?e+p^n + e+. (2) 

These rates depend on the known cross sections and the neutrino and antineutrino luminosities and mean energies. 
The cross section for Eq. (|^) is larger than that for Eq. (||) because of important weak magnetism and recoil 
corrections. 

The ratio of rates for Eqs. (|^) and (H) yields the initial electron fraction Ye 0, 

-1 

(3) 

Note as nuclei form. Ye will differ from Eq. (^). However, if this initial Ye is not significantly neutron rich, see below, 
a successful r-process is unlikely. In Eq. (|^) L is the luminosity, L the Ve luminosity and the mean Ve energy is, 

e^{El}/{E,J, (4) 

where the angle brackets indicate an average over the neutrino spectrum. We chose this definition for e because the 
lowest order cross section is proportional to {E'^ ) while the luminosity includes (E^^). We note, for a Boltzmann 
spectrum of temperature T, e = 4T = 4/3{E) while (E) = 3T. Likewise the mean Pg energy is. 



Ye 



^Q(e,e-)C(e,e-) 
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e={El)/{E,^). (5) 
The neutron proton mass difference A = M„ — Mp = 1.29 MeV or reaction Q value contributes the factor Q{e, e), 

The recoil and weak magnetism corrections to the cross section contribute the factor C(e, e), 

-'--)^^- 

Here M is the nucleon mass and, ao = 02 = {E^J(E^J/{EU^, describe spectral shapes. For simplicity 

we assume oq « 02 ~ 1.2 for both neutrinos and antineutrinos 0]. Although the factor C can increase Ye by 20%, it 
is neglected in many supernova simulations. 

Equation (^) involves observable quantities, the neutrino luminosities L, L and mean energies e and e. Therefore, 
one should be able to deduce Ye from observations of the next galactic supernova. This is a striking feature of r- 
process nucleosynthesis in neutrino driven winds. Perhaps the most important parameter Ye can be directly probed 
with observations. If observations suggest that Yg is not suitable for the r- process (see below), then the neutrino 
driven wind may be strongly disfavored as an r-process site. 

Supernova SN1987A already provides important data for Eq. (||). Jegerlehner, Neubig and Raffelt [|| place limits on 
the time averaged antineutrino temperature Tp^ from Kamikande and 1MB observations. We discuss time dependence 
below. If one neglects neutrino oscillations, Jegerlehner et al's 95% upper limit is, 

Tp^ < 4.6 MeV, (8) 

assuming a Maxwell-Boltzmann spectrum for which, 

e = ATp^ < 18.4 MeV. (9) 

If one assumes neutrino oscillations with parameters of the large mixing angle solution to the solar neutrino problem 
and if one assumes the muon or tau anti- neutrino temperatures are 1.7 times the electron anti-neutrino temperature 
then the limit decreases to § , 

Tp^ < 4.2 MeV, (10) 



e = 4Tp^ < 16.8 MeV. (11) 

Neutrino oscillations occur after the neutrinos pass through the wind and mix some hot i/^ or u^- into P,,. Therefore 
the original De spectrum must have been even colder so that the mixed spectrum could be consistent with observation. 
Note, Eq. ( pT| ) would be even lower but for matter effects as SN1987A neutrinos passed through the Earth. 

Originally, some groups assumed theoretical neutrino spectra and tried to use SN1987A data to set limits on 
oscillation parameters. Here, we assume that the oscillation parameters will be determined by independent (solar 
neutrino) experiments. Then, one should use the most accurate oscillation parameters to infer neutrino spectra from 
SN1987A data. We note that oscillations do not change our results very much. However, we include both oscillation 
and non oscillation cases for completeness. 

The ratio of anti-neutrino to neutrino luminosities is constrained by lepton number conservation. We assume, 

^<1-1- (12) 

Oscillations among active i/^, f^i, i't, flavors are not expected to greatly change this. However oscillations of electron 
neutrinos to sterile neutrinos (without charged current interactions) i^e Vs could increase L / L in the wind . 

Figure 1 shows Yg contours for different e and e assuming L = I.IL. The wind is neutron rich in the region to the 
upper left of the Y"e = 0.5 contour. The mean anti-neutrino energy e must lie within the shaded region, assuming no 
oscillations, or between the dot-dashed lines assuming oscillations. This sets a maximum value for e. In order for the 
wind to be neutron rich. 
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without oscillations and, 



e < 11.6 MeV, T^^ = | < 2.9 MeV, (13) 



e < 10.3 MeV, T^^ ^ ^ < 2.6 MeV, (14) 
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including oscillations. 

The requirement l^e < 0.5 is a reasonable minimum for the r-process. If one requires that the wind be significantly 
neutron rich, for example Yg < 0.4 than the limits become, 

e < 6.7 MeV, T^^ = ^ < 1.7 MeV, (15) 
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without oscillations and, 



e < 5.9 MeV, T^^ = ^ < 1.5 MeV, (16) 

including oscillations. The limits in Eqs. ( jJ^ through ( fJ^ j are significantly colder then most supernova simulations. 

Hoffman et al. discusses a range of physical conditions for the r-process. We consider two examples. The 
first scenario assumes a very short expansion time scale for the wind of order milliseconds ||l^,|l^ . This scenario can 
proceed with a relatively high ~ 0.48 because the number of seed nuclei formed is reduced by the short time scale. 
As a result, the ratio of free neutrons to seed nuclei can still be large enough to produce the heaviest elements. 

For this scenario, the limits in Eqs. ( [T3| ) and (|lj) are appropriate. However, the short expansion time scale may 
require a high neutrino luminosity that only occurs within a short time of core bounce. There may not be enough time 
for Ye to drop significantly in the protoneutron star. As a result, the opacity for Pe may not be drastically different 
from that for Ve and e may not be much smaller than e. Thus it may be difficult to satisfy Eqs. (|3|) or (p^ at short 
times. All realistic supernova simulations that we are aware of, for example [T5|-[l7t, do not satisfy Eqs. ( |13[ ) and ( |T^ ) 
at early times — say within 1/2 second of core bounce. 

A second scenario for the r-process involves a longer expansion time scale for the neutrino driven wind (of order 
a second). This can occur at later times in the supernova when the neutrino luminosities are lower. At later times 
there can be a large opacity difference between and so e can be significantly lower then e. 

However, the longer expansion time scale allows more seed nuclei to form. Therefore, one will need more free 
neutrons to have an acceptable ratio of neutrons to seeds. Thus Yg must be smaller, for example of order Ye < 0.4 



(assuming an entropy per baryon of order 130). The stringent limits in Eqs. ( p!5| ) and (16) may be appropriate for 
this longer expansion time scale r-process. These limits appear to be very hard to meet. Simulations tend to give 
much higher Ve temperatures. Again, we are not aware of any realistic simulation that satisfies them. 

We now discuss the time dependence of e. The cooling of the protoneutron star suggests that e should decrease 
with time. However, as the star becomes more neutron rich, the antineutrino sphere moves in, to higher densities 
and temperatures. This could cause e to increase with time. The SN1987A data suggest e decreases with time. The 
average electron energy for the six events detected at 1MB in the first five seconds is 36 MeV while the remaining 
two events at later times had an average energy of 20 MeV [Q. Likewise, the average electron energy of the eight 
events detected at Kamikande during the first five seconds is 17 MeV while the remaining three events had an average 
energy of 10 MeV Note, some of the difference in energy for the two detectors is explained by IMB's higher 

threshold. Clearly the data do not show a rapid increase of e at late time. However the statistics are poor and we 
have not reanalyzed the data to set quantitative limits on the time dependence. 

If the r-process occurs at early times, as might be necessary to have a very short expansion time scale, then the 
time dependence is not an issue. The bulk of the SN1987A events will directly constrain e. However, if the r-process 
occurs at late times, as is often assumed, then the time dependence could be more important. We emphasize that a 
very rapid increase of e with time, at late times, could evade our limits. However, there is no positive evidence from 
SN1987A for this rapid time dependence. 

Of course, one can evade our limits by assuming the r-process occurs in events with different neutrino spectra from 
those in SN1987A. This allows one to consider larger e. However, one may still have to explain a large ratio of e to 
e and why the events are different from SN1987A. Some alternative sites for the r-process include, prompt supernova 
explosions p^ , 'peculiar' supernovae jets |^ and colliding neutron stars [pT|-p3[. 

Equation (1^) includes weak magnetism and recoil corrections to the cross sections that are neglected in most 
supernova simulations. These are clearly important for Y^ in the wind. Therefore, they should be included in 
simulations. The reduction in the Ve cross section from weak magnetism may slightly raise the e or L predicted by 
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the simulations. However, e is constrained by SN1987A data while L/L is constrained by lepton number conservation. 
Simulations must still satisfy our limits on e independent of their inclusion of weak magnetism and recoil corrections. 

We use the Jegerlehner et al. results for e from SN1987A because of their simplicity. These results could depend 
somewhat on assumed spectral shapes, or other details. For example Janka and Hillebrandt |^ analyze SN1987A 
data assuming a Fermi Dirac spectrum with degeneracy -q, f{E) = [1 + e'^^/-^"''^]^^. They find 77 = is, somewhat 
weakly, favored with an upper bound of 2.5 and Tp^ < 4.5 MeV. For 77 = this is in good agreement with Eq. (||) 
while r] — 2.5 increases our limit in Eq. @ by 13% to e < 20.8 MeV. This, in turn, allows a 13% increase in the 
mean neutrino energy in Eq. ( p^ ) to e < 13.2 MeV. However, if the Ve spectrum has the same rj as the antineutrino 
spectrum, then the limit on the temperature in Eq. ([l^) is unchanged T^^ < 2.9 MeV. Thus, one can only increase 
the limit on T^^ by assuming a smaller 77 for v,, then for Pg. Note, we have not included small corrections from changes 
in ao and 02. The dependence of our bounds on spectral shape will be discussed further in future work. 

In conclusion: 

1) The electron fraction in the wind depends on observable neutrino luminosities and mean energies. Observations of 
the next galactic supernova should determine Ye- 

2) The mean antineutrino energy e is already limited by SN1987A data while the ratio of antineutrino to neutrino 
luminosities is constrained by lepton number conservation. 

3) To obtain a Ye suitable for an r-process in the neutrino driven wind, the mean electron neutrino energy could 
be significantly lower then that in present simulations. Alternatively, there could be a rapid increase of the mean 
antineutrino energy at late times. However, this is not seen in SN1987A data. Finally, there could be new neutrino 
physics such as oscillations to sterile neutrinos. 
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Mean Neutrino Energy (MeV) 

FIG. 1. Mean electron anti-neutrino energy e — {E'i^ ) / {Ei^^ } versus mean electron neutrino energy e — {E'i^ ) j (E^^ ) . Contours 
of constant electron fraction arc indicated for Ye values from 0.4 to 0.8. For a Boltzmann spectrum the electron neutrino 
temperature is 1/4 of e as indicated by the upper x-a^xis scale. The mean antineutrino energy e must be in the shaded region 
to be consistent with SN1987A data assuming no neutrino oscillations and between the dot dashed lines assuming oscillations 
with the large mixing angle MSW solar neutrino masses and mixing angles. 
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